Thermal Lifshitz Force between Atom and Conductor with Small Density of Carriers 
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A new theory describing the interaction between atoms and a conductor with small densities of 
current carriers is presented. The theory takes into account the penetration of the static component 
of the thermally fluctuating field in the conductor and generalizes the Lifshitz theory in the presence 
of a spatial dispersion. The equation obtained for the force describes the continuous crossover 
between the Lifshitz results for dielectrics and metals. 
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Introduction. The forces acting on atoms near the sur- 
face of a dielectric body was calculated at arbitrary dis- 
tances and temperatures by E. Lifshitz [1]. The Lifshitz 
theory is based on the calculation of the stress tensor 
for fluctuations of the electromagnetic field between two 
bodies. This tensor can be expressed in terms of the 
complex dielectric permittivities of the bodies at imag- 
inary values of frequency uj. The general equation for 
atom-surface forces was obtained as a limiting case of 
interaction between a dielectric body 1 with dielectric 
permittivity e (lu) and a body 2, in the limit when the 
second body is considered dilute from the electrodynamic 
point of view, i.e. assuming that dielectric function of the 
second body £2 (w) « 1 and expanding with respect to 
£2 — 1. In general, both zero-point and thermal fluctua- 
tions contribute to the force. In this letter we will discuss 
only the situation when the main contribution is due to 
the thermal fluctuations, as it occurs at hight enough 
temperature and large distance I between the atom and 
surface of the body. 

It was discovered by Lifshitz, that under the condition 



I » Afe^T = hc/ikBT) , 



(1) 



the interaction energy decreases according to the l/l^ 
law: 
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where a (0) and £ (0) are, correspondingly, the static val- 
ues of the electric polarizability of an atom and the di- 
electric permittivity of the body. 

This simple equation however, exhibits peculiar 
properties when applied to a conductor. Electrodynamic 
properties of a conductor at low frequencies can be de- 
scribed by a universal equation for the complex dielectric 
permittivity: 
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(3) 



where tr is the d.c. conductivity and e is the "bare" 
dielectric constant, which docs not take into account the 



contribution from current carriers. Because e (w) from 
([5]) tends to infinity as w ^ for any value of a, the 
potential tends to the universal limit 



^(0 = -^«(o) 



(4) 



which does not depend on e for an arbitrarily small con- 
ductivity - a seeming contradiction to common sense. 
However, the universality of ^ has a simple physical 
meaning. It expresses the fact that the static electric 
field does not penetrate into any conductor. Still, this 
statement, accepted in macroscopic electrodynamics, is 
only an approximate one. Actually, the field penetrates 
into conductors to a depth on the order of the well-known 
Debye radius Rd- For good conductors it typically is of 
the order of interatomic distances. However, if the den- 
sity of current carriers in the conductor is small, Rd can 
be large. We will see below that ([4]) is valid only at the 
condition 



Rd < I 



(5) 



(The problem of the temperature dependence of the 
forces is a subject of active discussion. See 0, 0] and 
references therein.) 

Atom-surface interaction. As a first step we will derive 
the Lifshitz expression for the atom-surface interaction 
directly from the Matsubara Green's function formalism 
of (see also Q). The basis of this theory is the equa- 
tion for the variation 5F of the free energy for a small 
change of the dielectric permittivity 5e : 



5F 



E' 

s=0 



Pf (C.;r,r)fc(i|G|,r)d3x , (6) 



where Vj^ is the Matsubara Green's function of the elec- 
tric field, = IsnksT /h and sign " ' " means that the 
s — term is taken with the coefficient 1/2. Notice now, 
that the presence of an atom at the point can be con- 
sidered as a small change of the dielectric permittivity 



5e{LLi, r) — A-Ka [uj) 5 [r — r^) 



(7) 
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Substitution of ([7]) into ([6]) gives the final expression for 
the energy of the atom-surface interaction 



Vil) 



^a(z|C.|)[2?zf (G;r,r')]r_,_ 



(8) 



s=0 



where 'Dfj^{C,s]Y,T') is the Matsubara Green's function 
of the electric field for a dielectric half-space. (The in- 
teraction can be analogously expressed in terms of the 
real-frequencies Green's functions, see Q.) Stress that 
(IH) is an exact equation. The problem is to calculate Vf^. 
in a proper approximation, in our case taking into ac- 
count spatial dispersion of dielectric. The function Vf^. 
satisfies the equation 



dA-5u^ + %e{i\Q,v)5, 



X. 



Ik 



(9) 



= iTTh^S.kS (r - r') 



In media with spatial dispersion, which will be considered 
below, e(i |^s| , r) must be understood as a linear opera- 
tor, actig on the variable r. Equation ([5]) allows one to 
obtain the general Lifshitz result for the atom-surface in- 
teraction in a much more direct way than in the original 
Lifshitz paper because it is much easier to solve (|9]) for a 
half-space than for the two-body geometry of [l| . We will 
see that by using ([5]) one can obtain the limiting equation 
(12) in a very simple way. The point is that relativistic re- 
tardation effects are not important under conditions ([T]). 
Indeed, the general condition of neglecting retardation, 
i.e. the condition of the quasi-stationarity of the field, is 



CO < c/l. 



(10) 



At the "Lifshitz distances" ([T|) only the s = term in 
([5]) is important. The problem then becomes purely 
static and one may calculate the Green's function, ne- 
glecting relativistic retardation effects. This problem will 
be solved in the next section. Notice that the above- 
mentioned peculiarity of the small cr limit exists only for 
the s = term. 

Longitudinal Green's function. As a first step we will 
separate the Green's function into its longitudinal and 
transverse parts, i.e. will present it in the form 



where 



(rot),,I?;i^ = Oand = 



ET 



(11) 



(12) 



The central point of the derivation is the statement 
that the smallness of retarding effects implies the inequal- 
ity 



1 2?/?.^ I « IvR.^^l 



(13) 



This means that only fluctuations of electrostatic nature 
are important. Indeed, a simple estimate from ([9|) gives 
2?^^ - HCs/cfv^^ < 2?^^. Particularly, the leading 
s — term in ([5]) is defined only by the longitudinal 
function Vff^^ (C;r,r'). 

To obtain an equation for the longitudinal part, let us 
apply the operator di to both sides of Neglecting the 
term of the order of QV^^ , we find 

d,[eDg^]^~A^hdi5{v-v') . (14) 

The first equation ([72]) can be satisfied identically if we 
introduce a scalar function Lp according to 

Dg^Dg^^hd^dl^. (15) 

(Such a function was used previously in [7,].) Substitution 
of (dni) into (HH) gives the equation for (p 



d, [ed,^] = -4^5 (r - r') 



(16) 



This is the equation for the potential of a unit charge 
placed at point r'. 

It is worth noting that one can derive these results 
by performing from the very beginning a gauge transfor- 
mation under the D-function of the 4-vector-potential. 
Above we used the gauge Q, where I?oo = ^^oi = 0. How- 
ever, in the problems where retardation effects are small, 
the "Coulomb" gauge, where Dpi = and diDik — 0, is 
more proper. In this gauge the main contribution gives 
-Dooi which actually is equal to — (ys (compare § 85, [5|). 

We can now express the interaction between an atom 
and body in the terms of the function Lp: 



V{l)^^a{Q) [a,9>(0;r,r')]. 



(17) 



This equation is valid for for any body. Therefore it de- 
scribes all effects of the body size and shape. The func- 
tion must satisfy usual electrostatic boundary condi- 
tions on the surface of the body. Then the boundary 
conditions for Vf^ will be also satisfied. If an atom in- 
teracts with a body, which occupies the half-space z < 0, 
the boundary conditions are 



(18) 



The electrostatic problem is solved in Problem 1 to 
7, H. The solution is: 



1 



1 



<P = 



1 



(z + z'f + ix-x'y 



-1/2 



(19) 



where r = {z, x}. 

The first term in p9|) does not depend on the presence 
of the dielectric and must be omitted. Differentiating 
the second term in order to calculate 2?;^^ (Cs;r,r') and 
going to the limit x ^ x', z z' — * we find after 
simple calculations 



Vi^ (0;r,r) = ft9i9,V(r,r') 



_e{0)-l 1 
'e{0) + 12P 



(20) 
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Equation ^ follows immediately after substitution of 
^ into (HZl). 

Notice that the results which we discussed are valid 
only in the state of full thermodynamic equilibrium. In 
particular, p7|) is valid only if the temperature of the 
body is equal to the temperature of the black-body ra- 
diation that falls on the body from the z > half-space. 
If these two temperatures are different, the interaction 
energy decays as I ^ oo according to a slow 1/P law (see 

Bad conductors. The calculation of the contribution 
of free carriers in the longitudinal Green's function de- 
mands a microscopic approach. We will assume that 
the gas of carriers is not degenerate, i. e. obeys the 
Boltzmann statistics and, for simplicity, that the carri- 
ers have unit charges ±e (e is the elementary charge). 
Then perturbation of the density of, for example, the 
positive carriers produced by the static potential (f is 
= n(+) (e"<='^/'=«^ - l) « -{n^+^eip)/{kBT). Ad- 
dition the corresponding charge density into (|16p gives, 
for 2: < 0, z' > 0, the following equation for ip: 



To solve the equation, let us expand (/^(x, z, z') into a 
Fourier integral with respect to x: 



ipk{z,z') = / (p(x, z,z')e 



'^■''d'^x 



(24) 



Let us consider first the z > domain. The Fourier 
transform of Eq. (fl6|) is 

(25) 



[d^^ - fc2](^k -inSiz - z') 



The presence of the (5-function means that ip\^ a,i z ^ z' 
must have the singularity of the type 27r|z — z'\. One 
easily finds the solution: 



-27r- 



-'k\z — z I 



+ Ae 



-kz 



(26) 



where A{z') must be defined from the boundary condi- 
tions. 

In the domain z' > 0, z < the equation for (/3k can be 
obtained by the Fourier transform of (PT|) : 



[dl - g2]^k = 



(27) 



[A 







where 



IkTf 



(21) 



(22) 



with 



(compare § 78, [9]). Here n = rS~^ -f rS'^^ is the total 
density of current carriers and e is the "bare" dielec- 
tric constant, without contributions from the carriers. 
Note that k = 1/Rd, where Rjj is the Debye radius. 
The Debye radius in good metals is of the order of inter- 
atomic distances. However, in "bad" conductors, where 
number of carriers is small, it can be comparable with 
the distance between the atom and surface. Equation 
(PT|) describes the screening of the electric field by the 
free carriers in the conductor. Indeed, according to this 
equation, the potential around a unit point charge in a 
uniform medium is e~'^^ /er. Notice that this equation 
can be interpreted in terms of spatial dispersion by in- 
troducing the longitudinal dielectric permittivity, which 
depends on the wave vector: e^(/c) — e[l + 1/ {kRn)'^]. 
Accounting for the spatial dispersion is a difficult prob- 
lem for restricted bodies. (See, for example [ll[ and a 
review [l2j and references therein. In these papers the 
effects of the spatial dispersion on interaction between 
conducting bodies are discussed.) It was easily solved in 
our static case due to the local connection between the 
electrostatic potential and the carrier density. 

Thus the function ip satisfies (ITB|) at z > and (1^ 
at z < 0. On the boundary z = we have now the 
"microscopic" boundary conditions 



(23) 



(28) 



The solution is 93k — Be"^^ . Substitution of the boundary 
conditions gives the equations for A and B: 



B = ~2Tre-''''' /k + A 
eqB 



-27re-'=^' - kA 



This gives 



^ , 2tt eq - k 

A{z') = -e 

^ ' k eq + k 



(29) 



(30) 



We can omit now the "free space" first term in (|26p and 
write the function (y9k as 



^k(z,.') = x^^''^'''''' 
k eq + k 



(31) 



Notice that ipi^ depends on z, z' only in the combination 
z + z'. This means that diO'i^p ~ [dl — d'^f, or for the 
Fourier components [9/9jV]k = ~2k'^ip\^. We can put now 
z = z' = Z,x = x'. Finally the potential of interaction 
(flT)) can be written as 

V{1) = ~fcBTa(0) / ^^^e-^'^^k^dk . (32) 
Jo eq + k 

If K which corresponds to the case of a "very 

bad" conductor, q k and we recover the result Q 
for an ideal dielectric, changing only e(0) —>■ e. In the 
opposite limit k ^ l^^ (|32|) is reduced to the "good" 
metal result (jlj. Equation ([32|l can be written as 



V{1) 



keT 
' 4/3 



(33) 
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where f = kI = {l/Ro) and Fq{^) is defined as 



H^dt . 



(34) 



From an experimental point of view the more interest- 
ing quantities are the force acting an atom —dzV and its 
derivative. They can be presented as 



-9iV^-^-^a{0)F,{0,d!V 



-a(0)F2(e) . 

(35) 

One can easily find expressions for these functions by dif- 
ferentiating Eq. (|3ip . According to the definition these 
functions 1 for ^ ^ oo and Fi (e — l)/(e + 1) 

for ^ — > 0. Functions Fi{^) are presented in Figure [T] for 
£ = 3.81 (fused Silica). One can notice the relatively slow 
convergence to the metal value at large 

Unfortunately the density of carriers in dielectrics de- 
pends very much on the technology of preparation of the 
samples and was not investigated systematically. For ex- 
ample, in [3] values of n in the interval between 10* and 
10^^ cm~'^ were used for fused silica. Then Ru is in the 
interval 2.3 x 10^ — 7.4 x lO^^^m. Prop erties of the vit- 
reous silica used in the experiments [l4llT5| are reviewed 
in an article [l6l |. This medium is an ionic conductor, 
the main carriers are ions Na"*". The density of Na im- 
purities is about UNa = 3 x 10^^ cm~'^, however it is 
difficult to estimate the number of ions which are effec- 
tive in mobility and screening. It is not presently clear 
if the presence of carriers has any influence on the in- 
terpretation of the results of the measurements 3, l3] ■ 
(This question was discussed in [l3].) Notice that the 
relaxation time of the charge distribution is very long 
in bodies with small a. This time can be calculated 
as Tc — e/(47rcr). At room temperature the resistivity 
p — 1/cr ^ 10^^ohm.cm= 1.1 x 10 ''s (see [17|). This gives 
Tc ^ 3.3 X 10^8= 917 hours. I believe, that at a such slow 
relaxation, the carriers mobility can hardly be important 
in any experiments. 

In conclusion, a theory for the interaction between an 
atom and a conductor due to the thermal fluctuations 
is developed. The theory takes into account the partial 
penetration of static electric fluctuations into the conduc- 
tor and is based on the Green's function technique in the 
presence of a spatial dispersion. A continuous crossover 
between an ideal dielectric and a good metal is investi- 
gated. 
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FIG. 1: EflFect of the field penetration in a silica sample on the 
atom-surface interaction. Functions -Fo(C) " solid line, Fi{£,) - 
dashed line and F2{^) - dotted line. 



[2: 

[s: 
[4: 

[s: 

[6: 

[7: 

[8 
[9 
[W 

[11 
[12 

[13 

[14 

[15 
[16 

[17 
[18 



[1] E. M. Lifshitz, Sov. Phys. JETP 2, 73 (1956). 



V. B. Bezerra et al., Phys. Rev. E 73, 028101 (2006). 
J. S. H0ye et al., J. Phys. A: Math. Gen. 39, 6031 (2006). 
I. E. Dzyaloshinskii and L. P. Pitaevskii, Sov. Phys. 
JETP 9, 1282 (1959). 

E. M. Lifshitz and L. P. Pitaevskii, Statistical Physics, 
Part 2, Pergamon Press, Oxford, 1984. 

C. Henkel, K. Joulain, J.-P. Mulet, and J.-J. Greffet, J. 
Opt. A Pure Appl. Opt. 4, S109 (2002). 

A. I. Volokitin and B. N. J. Persson, Phys. Rev. B 63, 
205404 (2001). 

L. D. Landau and E. M. Lifshitz, Electrodynamics of 
Continuous Media, Pergamon Press, Oxford, 1984. 
L. D. Landau and E. M. Lifshitz, Statistical Physics, Part 
1, Pergamon Press, Oxford, 1980. 

M. Antezza, L.P. Pitaevskii, and S. Stringari, Phys. Rev. 
Lett. 95, 113202 (2005). 

B. E. Sernelius, J. Phys. A: Math. Gen. 39, 6741 (2006). 
R Esquivel-Sirvent et. al., J. Phys. A: Math. Gen. 39, 
6323 (2006). 

A. C. Tien, S Backus, H. Kapteyn, M. Murnane, 
G. Mourou, Phys. Rev. Lett. 82, 3883 (1999). 

D. M. Harber, J. M. Obrecht, J. M. McGuirk, E. A. Cor- 
nell, Phys. Rev. A 72, 033610 (2005). 

J. M. Obrecht et al.,Phys. Rev. Lett. 98, 063201 (2007). 
Ch. M. Smith, in: Kirk-Othmer Encyclopedia of Chemi- 
cal Technology, , 5-th edition. Vol. 22, P. 407, John Wiley 
& Sons, Indianapolis, 2006. 

Fused Silica Code 7940 Brochure, Corning Glass Works, 
Corning, NY, 1992. 

G. L. Klimchitskaya and V. M. Mostepanenko, 
ArXiv:quant-ph/0801.1978. 



